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ABSTRACT 

Deep HST/NICMOS H (F160W) band observations of the z = 1.96 quasar 
pair Q 1634+267A,B reveal no signs of a lens galaxy to a la threshold of ~ 22.5 
mag. The minimum luminosity for a normal lens galaxy would be a 6L* galaxy at 
z ~ 0.5, which is 650 times greater than our detection threshold. Our observation 
constrains the infrared mass-to-light ratio of any putative, early-type, lens galaxy to 
{M/L)h ;^ 690/165 (1200/165) for = 0.1 (1.0) and Hq = 65/i65 km s"^ Mpc'^ We 
would expect to detect a galaxy somewhere in the field because of the very strong 
Mg II absorption lines aX z = 1.1262 in the Q 1634-1-267 A spectrum, but the HST 
H-band, I-band (F785LP) and V-band (F555W) images require that any associated 
galaxy be very under-luminous ^ 0.1L|^ (l.OLJ) if it lies within ^ 40/i^5^ (100/ig5^) kpc 
from Q 1634-F267 A and B. 

While the large image separation (3'.'85) and the lack of a lens galaxy strongly 
favor interpreting Q 1634-|-267A,B as a binary quasar system, the spectral similarity 
remains a puzzle. We estimate that at most 0.06% of randomly selected quasar pairs 
would have spectra as similar to each other as the spectra of Q 1634-1-267 A and B. 
Moreover, spectral similarities observed for the 14 known quasar pairs are significantly 
greater than would be expected for an equivalent sample of randomly selected field 
quasars. Depending on how strictly we define similarity, we estimate that only 0.01-3% 
of randomly drawn samples of 14 quasar pairs would have as many similar pairs as the 
observational sample. 
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1. INTRODUCTION 

The number of strong gravitational lenses has grown tremendously in the last five years, with 
well over 40 examples of strong gravitational lenses produced by galaxies (e.g. Keeton & Kochanek 
1996 and the CASTLES webpage When the separations are smaller than S'.'O, the lens galaxy 
has now been detected in all but one system: the small-separation, high-luminosity contrast lens 
Q 1208+1011. The lens galaxies have normal photometric properties and their mass-to-light 
ratios are typical of early-type galaxies embedded in dark matter halos (Kochanek et al. 1999b, 
Keeton et al. 1998). For the 18 systems with separations larger than S'.'O and nearly identical 
redshifts (see Table 1), a normal lens galaxy is detected in only four systems (RXJ 0911+0551, 
Q 0957+561, HE 1104-1805 and MG 2016+112). No lens galaxy has been detected in any of the 
remaining 14 systems, with typical lower bounds on the mass-to-light ratio of the lens about 10 
times higher than those measured for the normal lenses (Jackson et al. 1998). Whether the quasar 
pair^with nearly identical redshifts separated by S'.'O to lO'.'O are examples of "dark" gravitational 
lenses or binary quasars remains largely a mystery. Proving either scenario has very interesting 
ramifications on the dark mass distribution in the universe, or the nature of interacting quasars 
(Schneider 199S). There are about 2 such quasar pairs for every 1000 optically selected quasars 
(Hewett et al. 1998). 

Kochanek, Falco & Muhoz (1999a) used a comparison of the optical and radio properties of 
the quasar pairs to show that the majority of the systems must be binary quasars. First, five 
pairs (PKS 1145-071, HS 1216+50S2, Q 1S4S+2640, MGC 2214+S550 and FIRST J164S+S15) 
are "O^R" pairs in which only one of the quasars is a radio source and we can be certain that the 
system is a binary quasar. Second, even though most of the known lenses were radio selected and 
the angular completeness of the radio lens surveys is better than that of the optical lens surveys, 
there is no example of a "dark" radio lens (an "O^R^" pair). These two facts are inconsistent 
with the dark lens hypothesis. Quantitatively, Kochanek et al. (1999a) set a 2a (lex) upper limit 
of 22% (8%) on the fraction of the quasar pairs that could be gravitational lenses, the most likely 
scenario being that none is a gravitational lens. 

The high incidence and angular separations of the binary quasars can be quantitatively 



^http://cfa- www. harvard.edu/castles/ 

®The definition of quasar pairs also includes close physical systems which have been proven to be non-lenses, but 
exclude known lenses. 
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explained in terms of merger induced quasar activity (Kochanek et al. 1999a). The enhancement 
of the quasar-quasar correlation function by a factor of 10^ on the angular scales of the pairs is 
exactly as predicted by the locally observed enhancement by a factor of 10 in the quasar-galaxy 
correlation function on the same physical scales (Fisher et al. 1996, Yee &: Green 1987, French & 
Gunn 1983). The concentration of the pairs on scales smaller than lGf!0 (< 50 kpc) is a natural 
consequence of the need for a close passage and tidal interactions to trigger renewed quasar 
activity. The absence of binaries on scales smaller than 3'.'0 is a natural consequence of the rapid 
increase in the orbital decay rate due to dynamical friction as the binaries shrink. 

These limits on the existence of dark lenses are, however, statistical. While the maximum 
likelihood solution is that there arc no dark lenses, the statistical limit corresponds to allowing the 
existence of 3 dark lenses at the 2a limit. The existence of even 3 dark lenses would mean that 
the densities of dark lens halos and normal clusters are the same on these mass scales (Maoz et 
al. 1997, Kochanek et al. 1995, Wambsganss et al. 1995). Moreover, the binary hypothesis must 
somehow explain the remarkable similarity of the quasar pair spectra in the optical (e.g. Hawkins 
et al. 1997, Hewett et al. 1989, Steidel & Sargent 1991, Michalitsianos et al. 1997). 

Q 1634-|-267A,B remains one of the best candidates for exploring the issue of dark lenses 
and binary quasars. After the discovery by Sramek &; Weedman (1978) in a slitless spectroscopy 
survey, Djorgovski &; Spinrad (1984) and Turner ct al. (1988) found the optical spectra of the 
two quasars to be very similar, with an upper limit on the velocity difference of 150 km s^^, at 
the rcdshift oi z = 1.961. The flux ratio of the pair was 4.4 at R-band, although in images with 
seeing comparable to the image separation (Djorgovski &; Spinrad 1984). No lens galaxy has been 
detected down to limits of if 21.5 (Steidel Sz Sargent 1991) and R = 23.5 (Djorgovski &; Spinrad 
1984). 

Steidel &; Sargent (1991; hereafter SS91) used the Palomar 5-meter telescope to obtain very 
high signal-to-noise spectra of Q 1634-1-267 that illustrate the striking similarity of the spectra 
of components A and B, modulo a constant scaling factor of 3.28. They also found a number of 
absorption lines in the brighter component A, most notably the Mg II doublet AA2796,2803 at 
z = 1.1262, which perhaps hinted at the presence of a lens galaxy. When the emission line spectra 
of A and B were compared in detail, there were slight differences in the Ly a -I- N v. Si IV -|- 
O iv] , and C IV line profiles, and the line velocities were shifted by as much as 300-500 km s~^ 
(Small et al. 1997). The time delay created by the different travel times for the two rays can lead 
to spectral differences under the lens hypothesis, if the source quasar is variable. The temporal 
variations of quasar spectra have been examined by Filippenko (1989), Wisotzki et al. (1995), 
Impey et al. (1996), and Small et al. (1997). In particular. Small et al. (1997) demonstrated that 
the spectral differences observed in Q 1634-|-267A,B were typical of the differences seen in spectra 
of the same quasar taken after an interval of several years. 

With the advent of a second generation of instruments on the Hubble Space Telescope (HST), 
NICMOS has offered unique capabilities to discover and identify missing lens candidates by 
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observing them at rest optical wavelengths where they are bright, while providing unprecedented 
high angular resolution in the infrared necessary for studying these systems in detail. The 
CfA/Arizona Space Telescope Lens Survey (CASTLES) project is currently conducting a study 
of roughly 40 gravitational lens systems and candidates using optical/near infrared imaging with 
the HST. The selection sample consists of "simple" systems where the strong lensing is believed 
to be caused predominantly by a single galaxy. Some of the goals are: to place tighter constraints 
on the lens geometry; to study the mass and light distribution of the lens galaxy; to study 
the evolutionary history and environment of the lens galaxy; to better constrain cosmological 
parameters through lens modeling and statistics, and finally, in some cases to discover the missing 
lens galaxy. 

In Section 2 we present our observations of Q 1634+267A,B, the photometry, and the 
determination of the magnitude limit at the anticipated lens position. We also present new limits 
on the radio flux of Q 1634+267A,B and two other quasar pairs. Section 3 presents parameters 
derived for a hypothesized lens galaxy from a simple isothermal sphere model. We discuss the 
properties of the metal line absorber at z = L1262 in Section 4, and the statistical similarity of 
quasar spectra in Section 5. Conclusions follow in Section 6. We adopt Hq = 65 km s~^ Mpc~"^ 
throughout the paper and display key results for both Oq = 0.1 and Qq = 1. 

2. OBSERVATIONS AND ANALYSIS 
2.1. Observations 

We observed Q 1634+267 on 15 August 1997 with HST using NICMOS Camera 2 (NIC2) and 
the H (F160W) filter. We obtained four 690-second integrations in a 10.5-pixel dither pattern, 
for a total integration time of 2560 seconds. We reduced the data with our own software package 
"nicred" (McLeod 1997; Lchar ct al. 1999). Figure 1 shows the reduced, combined image of the 
field of Q 1634+267. Q 1634+267 A and B are approximately centered in the NIC2 field, and a 
star is visible at the edge of the ~ 19" x 19" field of view. 

2.2. Astrometry and Photometry 

The plate scales of the NICMOS cameras changed as the IR array underwent thermal 
expansion. The variation was monitored (Cox et al. 1998), and the NIC2 plate scales were 
0'.'0760926 and 0'.'0754090 per pixel in the X and Y directions, respectively for the measurement 
nearest the date of observation. The difference in the X and Y plate scales is due to a small tilt of 
the IR array relative to the focal plane. Wc adopt a zero point of A4h = 21.79 ± 0.02 for F160W 
for infinite apertures, based on a comparison of HST archival and ground-based observations of the 
standard star P330E (Persson et al. 1998). It agrees closely with the value of 21.83 mag obtained 
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by the Space Telescope Science Institute The foreground Galactic extinction in the direction of 
Q 1634+267 is only E{B - V) = 0.072 for Ry = 3.1 (Schlegel, Finkbeiner & Davis 1998); hence, 
we applied no extinction corrections to our photometric estimates. 

We fitted a point-spread function (PSF) to our combined image, to determine the relative 
coordinates and the magnitudes of Q 1634+267 A and B. Because A and B are well separated, 
we derived a PSF estimate from both components with the IRAF^ DAOPHOT package. Table 
2 shows the resulting relative coordinates and H-band magnitudes of A, B and the star within 
our field of view. The quasar images are separated by 3'.'852±0.003; the brightness ratio A/B is 
3.34 + 0.06. 

We obtained 10-minute snapshots of three quasar pairs, Q 1634+267, Q 1120+0195, and 
LBQS 1429-008, on 9 Dec 1998 with the VLA in the C array at A = 3.6 cm. We detected no radio 
flux to la limits of 0.028, 0.030, and 0.045 mJy/Beam, respectively. With these low flux limits, 
there is little hope that further radio observations can shed much light on these systems. 

2.3. Magnitude Limit at the Expected Lens Position 

Visual inspection of the images revealed no lens candidate near the expected position. To 
search for a faint lens that might be hidden by A or B, we fitted the PSF derived in the previous 
section to the quasar components. We then subtracted the fitted PSF from these components. 
The fainter one, B, was marginally broader than A. The difference might be caused by intrinsic 
PSF differences at different positions on the detector, and similar small PSF mismatches were 
seen in other NICMOS archival images containing two or more stars. Thus, we cannot conclude 
that the slight broadening is due to a faint galaxy underneath quasar B. Such a galaxy would have 
niH ^ 21.8 mag, which is too faint and at the wrong location to be able to produce lensed images 
with the observed splitting of ~ 3'.'8. 

To determine our detection limit for a lens galaxy, we generated circular galaxies with a de 
Vaucouleurs profile and an effective radius corresponding to an L* galaxy (~ Shg^ kpc) over the 
redshift range 0.4 ^ z ^ 1.8 and for = 0.1 and 1.0 cosmologies. We convolved the redshifted 
galaxies with the PSF derived from A and B, and added them to our NIC2 frame on a grid of 
positions within the frame. The resulting images were smoothed with a Gaussian kernel with cr = 3 
pixels, and inspected for signs of the added galaxies. We found a magnitude limit of mn = 22.5 
{la) for an early-type lens galaxy. For the sky background in our frame, the standard deviation of 
the noise level per pixel corresponds to a surface brightness of roughly 23.3 mag arcsec"^. 



^http: //www. stsci.edu/ftp/mstrument_news/NICMOS/NICMOS_phot/keywords. html 

*IRAF(Image Reduction and Analysis Facility) is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in Astronomy, Inc., under cooperative agreement 
with the National Science Foundation. 
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3. LOWER LIMIT ON A LENS M/L RATIO 

We compute the expected mass to light ratio for a putative lens galaxy for Q 1634+267 from 
our detection threshold. The anticipated mass of the lens can be derived from simple lens theory 
assuming it has a singular isothermal sphere (SIS) mass profile. The lens equation relates the 
angular coordinates ft {9) of the source (images), (5 = 6 — a{0), where a(6) is the deflection angle 
for an image at 9, with the lens galaxy at the center of coordinates. The SIS assumption implies: 
a{6) = (1/2)A0 = A7r{Dds/ Ds){cry'^ /c^), where Dg (Dds) is the distance from the observer (lens) 
to the source; is the observed angular separation of the lensed images, and is the velocity 
dispersion of the lens galaxy. Therefore, using the observed image separation one can determine 
the velocity dispersion for a lens galaxy as a function of redshift. We show in Figure 2a the 
velocity dispersion ay{z) using cosmologies Qq = 0.1 and 1, and Hq = 65 km s~^ Mpc~^. For an 
A/B magnification ratio of /i, the lens galaxy should lie A9{fj,/{fi + 1)} from image A along the 
A-B separation vector. 

Wc translate our magnitude limit oi H = 22.5 into an upper bound on the luminosity of the 
missing lens galaxy, placed at various rcdshifts. We follow the method implemented in Keeton et 
al. (1998) to derive iif-band luminosity of elliptical galaxies at various redshifts using Chariot &; 
Bruzual (1993) models. The models account for k-correction and passive evolution from formation 
redshift z = 5. Coupled with the masses required to produce lensing, determined via the lens 
equations, the luminosity limits yield the mass-to-light ratio (M/L) as a function of redshift, shown 
in Figure 2b for cosmologies of 0,q = 0.1 and 1. The most likely (comoving) distance to the lens is 
half the distance from the observer to the source, where the galaxy has the largest cross section for 
multiply imaging a source. For Q 1634-1-267 at redshift 1.961 the most probable lens location is 
therefore at redshift 0.45. An elliptical galaxy aX z = 0.5 would have M/L ^ 1000 in a low density 
universe, which is ^ 50 times higher than a normal lens galaxy (Keeton et al. 1998, Jackson et al. 
1998), and 50 times higher than a limit placed by Turner et al. (1988) in the optical. Therefore, 
if obscuration is responsible for hiding a lens galaxy with a normal mass to light ratio from view, 
then the dust must cause at least 4 magnitudes of extinction in the i7-band. If, however, we accept 
the Mg II absorption lines as evidence for a lens galaxy at z ~ 1.1 (SS91), then the detection limit 
of H > 22.5 corresponds to about 0.15L^, or M/Lh ^ 690 for an E/SO type galaxy. But if the 
lens galaxy type is an Sb, at a redshift of z = 1.1 it would have a M/Lh ^ 430. The minimum 
M/Lh is 400 for an Sb galaxy at a redshift of 2; = 1.3. 

If the lens is an early-type galaxy, we can predict its luminosity from the Faber-Jackson 
relation (1976) between velocity dispersion and luminosity, (Ty/cr* = (L/L*)^/^, where a* and L* 
arc the velocity dispersion and luminosity of an L* galaxy. Studies of lens statistics (e.g. Kochanek 
1996) show that the velocity dispersion estimated from the image separation should closely match 
the central velocity dispersion of the stars. Figure 2a shows the velocity dispersion estimated from 
the lens geometry as a function of the lens redshift. Kt z ^ 0.5 the predicted velocity dispersion 
is (7 !=a 350 km s"^, considerably higher than the characteristic velocity dispersion of an L* galaxy 
of cr* ~ 220km s~^. A galaxy at z pa 0.5 would be a L pa galaxy, far brighter than normal 
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early-type galaxies, where IL^ would correspond to H 17 A (16.9) mag for a passively evolving 
stellar population formed at ~ 5 using the Chariot & Bruzual (1993) models for Q = 0.1 (1.0). 
Such a bright galaxy would be at least 650 (1000) times above our detection limit. At the redshift 
of the Mg II absorber of z ~ 1.1, the required velocity dispersion rises to « 500 km s~^, at the 
highest end of the range observed for bright galaxies (see for example. Bender et al. 1996). It 
would have an extreme luminosity of 15L|^ if the scaling relation continued to hold. 

We conclude that the absolute lower limit on the mass-to-light ratio of the putative lens is 
M/Lh ~ 690 for an early type galaxy at a redshift z ^ 1.1, and M/Lh ^ 400 for an Sb type 
at a redshift of z 1.3. The limit is higher for any other redshift and for JIq > O-l- For the 
most probable lens redshift of z « 0.5 and the cosmology of a high-density universe, the limit 
is M/Lh > 1000. Depending on the redshift of the presumed galaxy, the observed H band 
corresponds roughly to rest R band {z 1.5) or / band {z « 0.4). Since {B — H) « 5.7 mag for 
ellipticals at redshift 0.5, and {B — H)q = 1.91 mag, the limit would be about 30 times higher for 
M/Lb than for M/Lh- If Q 1634-|-267A,B is the split image of a single quasar, the lens must be 
a galaxy with an unprecedentedly high mass-to-light ratio. 

In fact, as we discuss in §4, there are no galaxies above a detection limit of 0.1L|^ within 
40/i^5^ kpc and no galaxies above a detection limit of 1L| within lOOhg^ kpc of the quasars. 
Thus the depth of the imaging also rules out any normal cluster of galaxies centered near the 
expected lens position. Moreover, the limit on the optical M/L oi M/Lb ^ 10'^/i65 greatly exceeds 
even the usual estimates for groups and clusters of galaxies, where virial M/L k. 200-400/i (e.g., 
see Carlberg et al. 1996, Bahcall et al. 1995, David et al. 1995) is typical. 

4. THE NATURE OF THE Mg ii ABSORBER 

In Q 1634+267A,B, metal absorption lines are observed at z = 1.839 and z = 1.126 (Sargent 
& Steidel (1991), Djorgovski & Spinrad (1984), and Turner et al. (1988)). Mg II AA2796,2803 
absorption doublets with equivalent widths like the EWo(Mg II A2796) = 1.5lA doublet at 
z = 1.126 (SS91) are associated with chemically enriched material in the halos of visible, luminous 
galaxies along the line of sight (e.g. Steidel 1993, 1995; Stcidcl ct al. 1997). Because quasars 
show roughly one metal line absorption system per unit redshift down to a typical rest equivalent 
width limit of ~ 0.3A (e.g. Steidel 1993), the mere presence of a metal absorption system is not 
evidence for the presence of a lens galaxy. In a survey of galaxies selected by the presence of Mg II 
absorption, Steidel (1993) showed that the absorbers are associated with galaxies of Hubble types 
from E to Sb. The strength of the absorption is anti-correlated with the distance of the quasar 
from the nearest luminous galaxy, and it is correlated with the intrinsic brightness of the galaxy 
(Steidel 1995, 1993). The luminosities of the ~ 70 systems studied by Steidel (1995) and Steidel 
et al. (1997) range from 0.02L* to 4L* with 90% (61/70) having L > O.IL*. All but one (i.e. 99%) 
of the galaxies were found within a projected radius of 40/ig5^ kpc. 
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We searched both the NICMOS H image and existing HST WFPCl optical images (obtained 
under GTO-1116 by Westphal) for signs of a galaxy associated with the z = 1.126 absorption 
feature. Since the absorbing galaxy could have an exponential disk profile rather than the de 
Vaucouleurs profile considered in §2, we re-derived the detection limit for an exponential disk 
galaxy with a scale factor of 4 K^l kpc. The limit of 22.5 mag (Icj for fio = 0.1) was similar to 
that for a de Vaucouleurs model with a 3 K^l kpc effective radius. In the evolutionary models 
of Chariot & Bruzual (1993), a passively evolving L*, Sb-type, galaxy which formed at = 5 
would have ~ 19.7 for = 1-0, and ~ 20.0 for f^o = 0.1, at z = 1.126. Therefore, our 
detection limit implies that we could have detected a galaxy with L > 0.1L|^ anywhere within an 
impact parameter at AOhg^ kpc from Q 1634+267A (to the edge of the NICMOS array along the 
short direction). We derive a weaker detection limit of L > Lj from the F785LP (I-band) out to 
IOOHq^ kpc. Thus, any galaxy associated with the z = 1.126 Mg II absorption feature is either 
unusually faint or distant compared to those found in the Steidel (1993, 1995) surveys. 

Since the PSF of quasar B is broader than that of quasar A, it is unlikely that the null 
detection could be explained by a galaxy directly under Q 1634-I-267A, although we cannot 
exclude a very sub-luminous, compact dwarf with perfect alignment. Alternatively, a galaxy 
with a luminosity ^ lL|j lying outside our field of view might be responsible for the absorption. 
Distant (> 40 h^^ kpc), luminous (L > L*) galaxies are associated with absorption, but the Mg ii 
equivalent widths are less than 0.4A (Steidel 1995, Steidel et al. 1997), compared to the observed 
I.5A. For distant absorbers it is also odd that no equivalent Mg II absorption feature is seen in 
Q 1634-I-267B. In summary, the absorbing galaxy falls outside of the parameter space populated 
by galaxies found in Mg II absorber survey. The absorber must therefore be very under-luminous, 
^ O.IL^. 

5. SIMILARITY OF THE QUASAR SPECTRA 

In the absence of a lens galaxy, it is the remarkable spectral similarity of many of the quasar 
pairs which provides the impetus for the gravitational lens interpretation. The spectra of the 
Q 1634-|-267A,B pair show some of the greatest similarities of the pair population, and Small 
et al. (1997) have demonstrated that the observed differences are consistent with the temporal 
variations of quasar spectra over epochs separated by months to years. No study, however, has 
examined whether the spectra of the close pairs show a degree of spectral similarity that is 
statistically greater than that of randomly selected quasars. An affirmative answer increases the 
credibility of the lens hypothesis, because the lens hypothesis provides a natural explanation for 
the spectral similarities. While it would not rule out the binary hypothesis, it leaves a serious 
puzzle as to why the nuclear regions of two separate galaxies have such highly correlated emission 
properties. We start by defining how we will measure spectral similarities, and then estimate the 
likelihood of finding a single quasar pair with spectra as similar as Q 1634+267A,B. In isolation, 
the significance is hard to interpret because we picked Q 1634+267A,B for study due to its spectral 
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similarities. Thus, in our final step we estimate the Ukehhood of finding a population of quasar 
pairs with the observed spectral similarities. 

5.1. Sample Definition 

We will compare spectra using the observed differences between the spectral slope a and 
the rest equivalent widths (EWq) of the Ly a, C IV (AA1548, 1550) and C iii] (A1909) emission 
lines. We selected these features because they are commonly measured in statistical surveys of 
quasar spectra. For Q 1634+267 we computed the indices from the SS91 spectra. These calibrated 
spectra were taken at the Palomar 5-meter, using a double spectrograph that covered the blue 
(3150-4750A) and the red (4600-7000A) wavelengths. We used either the Large Bright Quasar 
Survey (LBQS) or the pair population as a whole for a comparison sample. 

The LBQS consists of 1055 optically selected quasars (Foltz et al. 1987, Foltz et al. 1989, 
Hewett et al. 1991). We use the measurements of the spectral indices for the LBQS quasars from 
Francis (1996), Hewett et al. (1995), Francis (1993a), and Francis et al. (1991). The binned 
distributions for the LBQS quasars were fitted with analytic functions (composites of Gaussians 
and Lorentzians) to provide smooth probability distributions. Some of the spectral features are 
correlated (Francis et al. 1992; Francis 1993a), and we will conduct tests to examine the effects of 
these correlations on the probability calculations. 

Unfortunately we lack a complete spectral database for the pairs, although the uncanny 
similarity of Q 1634+267A,B is not an isolated case. There are 14 quasar pairs with angular 
separations 3" ^ d ^ 10" and similar redshifts (see Table 1). We will frequently refer to this 
as our 14 pair sample. Of these 14, we can reject 7 as lenses. In the cases of PKS 1145-071, 
HS 1216+5032, Q 1343+2640, J 1643+3156 and MGC 2214+3550, only one of the quasars is a 
radio source; in the case of Q 0151+048 the small emission line redshift difference is confirmed by 
an absorption feature in the background quasar; in LBQS 2153-2056, the spectrum of quasar B 
shows a strong absorption blueward of the C IV emission line. Four of these binaries (QJ 0240-343, 
PKS 1141-071, HS 1216+5032, Q 0151+048) show fairly similar spectra, although they ah have at 
least one significant spectral difference. Of the 7 ambiguous pairs, five (Q 2138-431 in Hawkins et 
al. 1997, Q 1429-008 in Hewett et al. 1989, Q 1634+267 and Q 2345+007 in SS91, Q 1120+0195 
in Michalitsianos et al. 1997) show remarkably similar optical spectra. For the pair sample we 
restrict our comparison to the C IV equivalent width and the power law slope a because the other 
two indices are not consistently observed or are impossible to measure without the original data. 

5.1.1. Line and Continuum Slope Measurements 

To determine the spectral slope, we fit a power law (defined as F,y oc i/") excluding regions 
where there are emission or absorption lines. The Fe II in the UV spectra of quasars introduce 
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substantial uncertainty in the measurement of continuum slopes. The power law a is also uncertain 

in the blue spectrum from contamination due to the extended wings of the C IV, N v, and Ly a 
emission lines, as well as the Ly ol forest absorption. Because of these complications, the width of 
the spectral slope distribution remains controversial. Some (e.g. Elvis et al. 1994; Webster et al. 
1995) find large dispersions, while others find little or no dispersion in a (Francis 1996; Sanders 
et al. 1989). In our analysis we will adopt the narrow a distribution from Francis (1996), as 
derived from 30 LBQS quasars. The Francis (1996) distribution has a root-mean-squared (RMS) 
width of cTq = 0.3, compared to a a = 0.5 for the broader distributions of Francis et al. (1992). In 
Q 1634-1-267 the mean slopes are a = —0.17 in the blue (observed 3150-4734A) and a = —0.12 in 
the red (observed 4580-6985A) spectra respectively. In both the blue and the red, the power law 
indices between A and B differ by Act = 0.03 ± 0.03. 

When comparing equivalent widths we must be careful to use identical measurement methods 
because of the difficulties introduced by line blending and continuum definition. Francis (1993a) 
showed that C IV equivalent widths can vary by a factor of two depending on the measurement 
method, and thus broaden the distribution of equivalent widths by the same factor. For 
Q 1634-|-267A,B we avoid these problems by using the same definitions followed by the LBQS 
(Francis 1993a). For this method, continuum windows are defined on both sides of a given line. 
Then the flux is summed above a line adjoining the continuum windows, over a region known as 
the line window. The continuum windows are chosen to exclude the flux contributed by nearby 
broad lines, so the line fluxes should be considered a lower limit. The difference in rest equivalent 
widths (EWo) between A and B arc SA for C IV, 8A for Ly a and 3A for C III]. Measuring 
EWo(C iv) requires connecting the red and blue halves of the spectra because the blue side 
spectrum terminates on the red wing of the C IV emission line, causing the continuum windows to 
fall on different spectrographs. We calculate a scaling factor of 1.22 using the regions of overlap. 
But there is some evidence that the same scaling factor might not apply to quasar B because of 
large fluctuations in the data over that region. Therefore the EWo(C iv) differences of 3A should 
be considered an upper limit. Indeed, restricting ourselves to the blue regions alone and fitting 
Gaussians to the emission lines we find quasars A and B to differ at most by AEWo(C iv) = lA. 
The two methods yield the same flux. For completeness, and because the numbers have not been 
published elsewhere, we also measure other line fluxes and summarize them in Table 2, although 
some we do not use in our analysis. All of the line fluxes are measured with the technique 
described in Francis (1993a) except for N v. Because the blending with Ly a, the measurement 
for N V lines is done by fitting Gaussians to Ly a simultaneously with N v where we assumed 
that both lines are similar in shape. 

For the other quasar pairs we estimated bounds on the C IV equivalent width and the spectral 
index a based on the published spectra. For the 5 most similar pairs, we estimate that they have 
spectral differences Aa ^0.1 and AEWo(C iv) <^ 6A (see Table 3). For the two pairs whose 
spectra have the worst noise (Q 1120-^0195 (Michalitsianos et al. 1997) and Q 1429-008 (Hewett 
et al. 1989)) the estimates are rough, while they are more accurate for the other three pairs. Since 
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the estimates of the equivalent widths may be imprecise, we will also quote results by doubling 
our bound on the differences to AEWo(C iv) ^ 12 A. We summarize our parameter estimates for 
the pairs in Table 3. 

5.1.2. Joint Probabilities and Correlations among Spectral Features 

To compute the probability that two random quasars have given differences in C III], Ly a, 
C IV equivalent widths and spectral slope a, we randomly draw the properties of two quasars at a 
time from the compiled number distributions. If the spectral features are completely uncorrelated, 
the likelihood that two quasars are similar to each other, i.e. have differences of less than 
AEWo(C iv), AEWo(C lll]), etc., is simply the joint product of each individual probability. 
However, if there are correlations, they could make the alarming similarities of the quasar pairs 
a trivial selection effect. Therefore we first explore the extent to which spectral features are 
correlated. 

Through a principal component analysis, Francis et al. (1992) showed that the equivalent 
widths of quasar emission lines are correlated with the power law slope a. Francis (1993b) further 
showed that a is correlated with redshift. Their study was based on correlations with a spectral 
slope distribution in Francis et al. (1992) which has a large RMS dispersion of about ctq = 0.5. The 
strongest correlation, Al lll+C III] with a, has a large scatter; the correlations of this and other 
lines are very weak. Subsequently, Francis (1996), using more accurate photometry, found that 
the spectral slope distribution actually has a much narrower dispersion of aa = 0.3. But the newer 
study does not affect the dispersions of the line strengths, which are more accurate as long as 
they are measured consistently. Because the scatter of the line strengths remains large for a much 
reduced range in spectral slopes, the correlations are weaker in light of the narrow distribution 
compared to Francis et al. (1992). Furthermore, plotting spectral slope vs. redshift from Francis 
(1996) shows that there is no correlation between the two parameters. In the analysis below, we 
will use the narrow power law distribution found in Francis (1996). Therefore, the correlations 
between line emission, continuum, slopes and redshifts can be safely ignored for our purposes. 

Nonetheless, to confirm that correlations are unimportant even under extreme scenarios we 
conduct the following Monte Carlo experiment to choose quasars from the number distributions 
found in Francis et al. (1992) and Francis (1993a). We select from the broader spectral slope 
distribution where correlations were originally found. The correlation between emission line 
strengths and a, and a with redshift z, means that the process of drawing quasars must proceed 
in a sequence to account for correlations at each step: first we draw quasars out of the LBQS 
redshift distribution, then the spectral slope distribution, followed at last by the emission line 
distribution. To account for the correlation of a with redshift we do the following. Having picked 
a redshift, and given the power law distribution a, we offset the entire distribution by an amount 
da = 0.79 X d log{l + z) (See Figure 2 of Francis 1993b). We can then randomly select a value 
of a from this new distribution. There is, however, one subtlety. The distribution for a contains 
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no redshift information because it is summed over all redshifts. But because of the correlation, if 

quasars at any given redshift have a small dispersion in their spectral slopes, as one accumulates 
the distribution over broader redshift ranges, the net distribution would also gradually become 
wider. We account for this effect by arbitrarily shrinking the dispersion and then drawing from 
this narrower spectral slope distribution instead of the original, broader one. We verified that 
decreasing the dispersion by 25%, more than one would expect, would not significantly alter our 
discussions below. Once a value for a has been selected, one can then proceed in the same fashion 
to select from the number distribution of C III], C IV, or Ly a, accounting for their roughly linear 
correlations with a by: EWune = ma + b where the correlation coefficients m and b are given 
in Francis et al. (1992). Again, we verified that any correlations would only begin to affect our 
conclusions if the number distributions were all narrowed by more than 25%. 

There may also be weak correlations of quasar spectral properties with luminosity. The pair 
members share a common redshift and have similar luminosities, so we would under-estimate 
the likelihood of the members possessing similar spectra if the distributions of spectral indices 
are correlated with either variable. We saw that the spectral index correlation with redshift 
has little or no role in our analysis. The Baldwin (Baldwin 1977) effect, a correlation between 
luminosity and C IV equivalent width, has been studied extensively (e.g. more recently see 
Kinney et al. 1990, Sargent et al. 1989, Baldwin et al. 1989). Their results show that at a given 
luminosity the equivalent widths of the C IV emission line have a defined range in strengths, with 
a scatter of about 20A. The strength of the C IV lines decreases with luminosity and the relation 
holds over 7 orders of magnitude in quasar luminosity. The parameterization of the correlation: 
EW{C iv) oc 7.1450 has /? ranging from (3 = -0.07 (Francis et al. 1992) to /3 = -0.45 (Smetanka et 
al. 1991). LBQS quasars, being a magnitude limited sample, have only a small spread in absolute 
luminosity, —28 ^ Mb ^ —26 (Francis et al. 1992). Since at any luminosity there is a considerable 
scatter in C IV and vice versa (see Osmer & Shields 1988 for summary), we do not expect the 
luminosity correlations to affect our analysis. Nonetheless, we tested for any effects by deriving 
the distribution of AEWo(C iv) both with and without limits on the luminosity difference AL of 
the quasars making up the simulated pairs using the Bj magnitudes and EWo(C iv) values for 
the 439 LBQS survey quasars from Francis (1993). We found no significant differences between 
the distributions without a luminosity restriction and those where the quasars were required to 
have absolute luminosities differing by less than a factor of 5 (1.7 mag) or a factor of 2.5 (1 mag). 

Because the parameter correlations are very weak, in the discussions that follow, we make the 
assumption that no two spectral features are correlated. 

5.2. The Spectral Similarity of Q 1634+267A,B 

As discussed previously, the spectra of Q 1634+267 A and B have equivalent width differences 
A < 3A for C IV; 8A for Ly a; and 3A for C III], as well as 0.03 for a. If we select randomly from 
the distributions of these features in the LBQS survey, the probability of finding two quasars with 
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smaller differences are 11%, 16%, 30%, and 12% respectively. Assuming no correlations between 
the features, the joint probability of the differences being as small as observed is 0.06%. 

With such a low joint probability, it would be reasonable to find several quasars with such 
similar spectra in the total quasar population (~ 10^ quasars), but highly unlikely to find any in 
the tiny sub-population of quasar pairs (14 objects). The strength of this conclusion depends only 
weakly on the details of the comparison, and the joint probability only rises to 1% if we confine 
the comparison to the spectral index and C iv equivalent width differences. Adding additional 
spectral features to the estimate would further reduce the likelihood, but it becomes increasingly 
important to fully understand the spectral correlations of quasars. Moreover, calculations of 
a posteriori probabilities can generally yield arbitrarily low likelihoods by over-specifying the 
problem. Nonetheless, the similarity of the Q 1634+267 A and B spectra extends across the entire 
rest-UV spectra of the pair. 



5.3. The Spectral Similarity of the Pair Population 

Because we pre-selected Q 1634+267A,B for study because of its spectral similarities, we 
cannot interpret it in isolation from the rest of the pair population. Moreover, the uncanny 
similarities of Q 1634+267A,B are not an isolated example among the pairs. To fully understand 
the significance of the spectral similarities of the pairs, we must consider the entire pair population. 

Under the assumption that the quasars making up the pair sample are drawn from the same 
parent population as the quasars in the LBQS, we can produce mock samples of quasar pairs by 
randomly selecting 14 pairs of quasars from the LBQS. We can then compare the distribution of 
the spectral slope and C IV equivalent width differences, Aa and AEWo(C iv), to those observed 
for the pairs. We will estimate the probability of finding at least as many similar pairs as observed 
in the real sample. 

We generated IC^ mock pair catalogs to determine the fraction of catalogs in which N^j^ 
pairs have spectral differences smaller than the specified limits on Aa and AEWo(C iv), with the 
results illustrated in Figure 3. The top panel shows the limits on spectral similarity (Aa ^0.1 
and AEWo(C iv)^ 6A) which characterize the differences observed for the 5 pairs with similar 
spectra (Q 2138-431, LBQS 1429-008, Q 1634+267, Q 1120+0195 and Q 2345+007). The 
chance of drawing 5 or more pairs with such similar spectra in the sample of 14 pairs is about 
0.2%. If we were to relax the constraints to Aa ^ 0.1 and AEWo(C iv)^ 12A to account for 
possible systematic errors in measuring AEWo(C iv), the probability for finding 5 or more such 
pairs grows only to 3%. We conclude that the likelihood for finding 5 pairs with such similar 
spectra in a sample of 14 pairs is low. If we significantly broaden the criterion for similarity 
(Aa ^ 0.2 and AEWo(C iv)^ 20A), the hkelihood of finding five or more similar pairs rises to 
60%. However, with these weaker criteria there are now seven observed pairs meeting the criterion 
for similarity (we should add QJ 0240-343 and LBQS 2153-2056), and the likelihood of finding 
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seven or more similar pairs is then smaller, at 20%. We regard this third case as a significant 
under-estimate of the similarities. The low probability for finding similar quasars randomly, 
using two parameters, may only be an upper bound. If our pairs sample were partly made up of 
gravitational lenses, microlensing and time delay might cause some spectral differences in those 
systems. It is impossible to quantify such a conjecture with our pairs sample. 

The probability for finding similar quasars is significantly lower if we use a broad continuum 
slope distribution (e.g. Francis 1992) rather than the narrow distribution of Francis (1996) used for 
the estimates in Figure 3. For the broad estimate of the a distribution, the likelihoods for drawing 
at least 5 similar pairs under the same three sets of constraints on Aa and AEWo(C iv) are 0.01, 
0.2% and 1.2% instead of 0.2%, 3.0% and 20% (see Figure 4). Thus, the debate about the width 
of the continuum slope distribution does not alter the conclusion that the spectral similarities of 
the pairs are unlikely to occur by chance unless the distribution is even narrower than that found 
by Francis (1996). The distribution of C IV equivalent widths, as long as measured consistently, 
does not suffer from the systematic problems of the spectral slope distribution. 

6. CONCLUSIONS 

We again confirm that there is no sign of a lensing galaxy or a surrounding group or cluster of 
galaxies near the Q 1634+267A,B quasar pair. Our limit on the mass-to-light ratio is ^ 690/i65 in 

the H band, tighter than previous estimates. We also lack a candidate galaxy within lOOh'^r^ kpc 
for producing the z = 1.126 Mg II absorption feature in the spectrum of A. Since the existence of 
even one large separation lens produced by a "dark" mass distribution has profound implications, 
confirming or disproving the lens hypothesis is highly desirable. It appears from the data in 
the literature that the flux ratio of the quasars is significantly time variable, changing from 4.4 
(Djorgovski & Spinrad 1984) to 2.83 (Turner et al. 1988), to 3.28 (SS91) and 3.34 (this work). 
With its large angular separation, regular photometric monitoring of the system to search for a 
time delay would be simple on modest size, ground-based telescopes. Although we found that the 
pairs Q 1634+267, Q 1120+0195, and LBQS 1429-008, are at best weak radio sources, with upper 
limits on their fluxes of 0.28, 0.030, and 0.045 mJy/Beam at A = 3.6 cm, deeper observations 
may still detect the sources or demonstrate that both members are at the extremes of the radio 
to optical flux ratio distribution. Deep X-ray observations would also be useful, both to measure 
the flux ratio of the quasars and to search for hot X-ray emitting gas associated with any lensing 
potential. Deep IR imaging should detect the quasar host galaxies, whose shapes must show 
tangential, arc distortions under the lens hypothesis. 

The spectral similarities of some of the quasar pairs has always been the major impetus for 
the "dark" lens interpretation of the systems. To date, these arguments have been based on the 
similarities of individual pairs rather than the overall pair population. Using the distribution of 
continuum slopes and C IV equivalent widths, we demonstrated that it was improbable to find as 
many similar pairs as observed in the existing sample of 14 pairs. The probability estimates range 



-15- 



from 0.01% to 3%, depending on how strictly similarity is defined and depending on whether a 
broad or narrow distribution was used for the continuum slope distribution. A better quantitative 
analysis needs a uniform spectral survey of the pairs. 

The effects of binary quasar interaction on spectral similarity have been surmised, for 
example, by Schneider (1993). Kochanek et al. (1999) demonstrated that most of the quasar 
pairs are binary quasars rather than "dark lenses." The binary quasar hypothesis is a successful 
quantitative interpretation of both the number of quasar pairs and the distribution of their 
separations. If they are all binary quasars, then there must be a physical explanation for the 
spectral similarities. Binary quasars are examples of quasar activity triggered by mergers, which 
represent a small fraction of all quasar activity (about 5% in simple models, see Kochanek et 
al. 1999a). We present two hypotheses for how binary quasars can have such similar spectra. 
First, if the physical properties of quasars triggered by major mergers (e.g. mass, accretion rate, 
surrounding gas density) were confined to a limited range of values compared to the overall quasar 
population, then the observable properties of the binaries would be more similar than those of 
randomly selected quasars. This hypothesis has a definite, testable prediction - the spectra of 
the different pairs should be more similar to each other than for randomly selected quasars. It 
would also be expected that the quasar host galaxies have properties commensurate with being 
the products of major mergers - which can be tested by deep imaging. Second, if quasar spectra 
have a significant dependence on the age of the activity, then the binary spectra are more similar 
than for random quasars because the merger provides a common triggering event for the activity 
in the two systems. Here we would not expect the spectra of the different pairs to show any more 
similarity than randomly selected quasars. In either case, the binary quasar population becomes 
a very important probe of the physics of quasars because it provides evidence for correlations 
between events on large scales and the detailed operation of the central engine. 
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Table 1. Quasar Pairs With Separations of 3" to 10' 



Name 




A6» 


R 


fo 




F20 
mjy 


|A^| 
km s ^ 


Lens? 


Type 


Ref 




Q 0151+048t 


1.91 


3'.'3 


28 


27.5 




< 1 


520 ± 160 


No 


02 


1 


PKS 1145-071 


1.35 


4" 2 


36 


2.1 


> 500 


740 


200 ± 110 


No 




2 


HS 1216+5032 


1.45 


9"1 


78 


5.2 


> 28 


3.9 


260 ± 1000 


No 


O'^R 


3 


Q 1343+2640 


2.03 


9'.' 5 


78 


1.1 


> 57 


8.6 


120 ± 890 


No 


O'^R 


4 


J 1643+3156 


0.59 


2'.'3 


17 


1.3 


> 40 


120 


80 + 10 


No 


O'^R 


5 


LBQS 2153-2056 


1.85 


7'.'8 


64 


14.5 




< 1 


1100 ± 1500 


No 


02 


6 


MGC 2214+3550 


0.88 


3'.'0 


26 


1.6 


> 42 


246 


148 ± 420 


No 


O^R 


7 


MG 0023+171 


0.95 


4" 8 


40 


3.0 


~ 10 


186 


292 ± 260 


?- 


O^R'^ 


8 


Q 1120+0195tt 


1.46 


6"5 


56 


173.8 




< 1 


628 ± 120 


?_ 


02 


9 


LBQS 1429-008 


2.08 


5{'l 


42 


nA 




< 1 


260 ± 300 


?_ 


02 


10 


QJ 0240-343 


1.41 


e'.'i 


52 


2.1 




< 1 


250 ± 180 


? 


02 


11 


Q 1634+267 


1.96 


3'.'8 


32 


4.4 




< 1 


33 + 86 


? 


02 


12 


Q 2138-431 


1.64 


4'.' 5 


38 


3.0 






+ 115 


? 


02 


13 


Q 2345+007 


2.15 


7'.'3 


58 


4.0 


> 9 


0.035 


476 ± 500 


? 


02 


14 


RXJ 0911.4+0551 


2.80 


.3"1 


24 


1.9 




< 1 


158 ± 1000 


Yes 


02 


15 


Q 0957+561 


1.41 


e'.'i 


52 


1.4 


1.3 


552 


200 ± 15 


Yes 




16 


HE 1104-1805 


2.32 


3'.'1 


24 


4.8 




< 1 


300 ± 90 


Yes 


02 


17 


MG 2016+112 


3.27 


3'.'6 


26 


1.7 


~ 1 


191 


40 ± 100 


Yes 


02 R2 


18 



Note. — Zs is the source redshift, AS is the angular separation, R is the projected separation at the source 
redshift for Qq = 1 and Hq = 50/i6o km Mpc~^, /o and fji are the optical and radio flux ratios or their 
lower hmits, Fx is the radio flux of the brighter image at wavelength A (in cm), and | Ad| is the velocity difference 
between the quasars. The entries in the Lens? column are: "Yes" if a normal lens (galaxy, group, or cluster) is 
seen in the correct position to produce the observed system, there is no significant velocity difference, and the 
radio and optical data are consistent with the lens hypothesis; "No" if we see no lens and either the radio emission 
or the emission line velocity difference, confirmed by an absorption line velocity difference, are inconsistent with 
the lens hypothesis; and, "?" if we see no lensing object but have no objective criterion to decide whether or not 
the object is lensed. If there is some evidence that the system is actually a binary, we used the label "?-" . Type 
denotes the optical/radio classification of the pair. Note that MG 2016+112 is really a triple system, not a pair. 
tQ 0151+048 is also named PHL 1222 and UM 144. 
ttQ 1120+019 is also named UM 425. 

References: (1) Meylan et al. 1990, (2) Djorgovski et al. 1987, (3) Hagen et al. 1996, (4) Crampton et al. 1988, 
(5) Brotherton et al. 1999, (6) Hewett et al. 1997, (7) Muiioz et al. 1997, (8) Hewitt et al. 1987, (9) Meylan & 
Djorgovski 1989, (10) Hewett et al. 1989, (11) Tinney 1995, (12) Djorgovski & Spinrad 1984, (13) Hawkins et al. 
1997, (14) Weedman et al. 1982, (15) Bade et al. 1997, (16) Walsh et al. 1979, (17) Wisotzki et al. 1993, (18) 
Lawrence et al. 1984. 
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Table 2. Observed Properties of Q 1634+267 A,B 



Parameter 


Value 




Q'/oco _|_ n nriQ 


A RA 


-0'/673 


A DEC 


3'/793 


H-magnitude of component A 


18.26 ±0.01 


H-magnitude of component B 


19.56 ±0.01 


H-magnitude of star 


19.67 ±0.01 


Flux ratio (A/B) 


3.34 ±0.06 


Power law slope 


-0.13 


EWo(Ly a) 


47.95A 


EWo(N V A1240) 


24.18A 


EWo(Si IV A1400) 


11.83A 


EWo(C IV AA1548, 1550) 


37.44A 


EWo(C III] A1909) 


16.91A 


^Power law slope defined as Fj, 


OC I' 



Table 3. QSOs A,B Parameter Difference Estimates 





Ao 


AEWo(C IV) A 


Lens? 


References 


Q 1120+0195 


0.1 


6 


? 


Michalitsianos et al. 1997 


LBQS 1429-0053 


0.03 


6 


? 


Hewett et aL 1989 


Q 1634+267 


0.03 


3 


? 


Steidel & Sargent 1991 


Q 2138-431 


0.02 


6 


? 


Hawkins et al. 1997 


Q 2345+007 


0.07 


3 


7 


Steidel & Sargent 1991 


LBQS 2153-2056 


0.05 


14 


No 


Hewett et al. 1998 


QJ 240-343 


0.2 a 


20 


No 


Tinney 1995 



Note. Estimated differences between the spectra of QSO pairs A and B from 
published spectra in literature. 

*Prom Tinney 1995. 
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Fig. 1. — Image of Q 1634+267 A and B. Horizontal bar shown is 2" in length, which corresponds 
to roughly 16 z /igg^kpc, for an f2 = 0.1 open cosmology. 
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Fig. 2. — (a) Anticipated velocity dispersion of a putative lens galaxy and (b) M/Lh ratio as 
function of redshift for a putative lens galaxy, assuming properties of an early-type galaxy. 
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Fig. 3. — The fraction of trials in which Ngim pairs from a sample of 14 pairs randomly drawn 
from the LBQS meet the stated conditions on Aa and AEWo(C iv). We assume the narrow 
spectral slope distribution from Francis (1996). The arrow marks the number of pairs observed to 
meet the similarity condition, (a) Representative spectral differences actually observed for the five 
lens candidates. Probabilities for (b) and (c) have substantially weaker C iv and spectral slope 
constraints, as indicated in each box. 



-24- 




Fig. 4. — The fraction of trials in which Nsim pairs from a sample of 14 pairs randomly drawn from 
the LBQS meet the stated conditions on Aa and AEWo(C iv). Here we assume the broad spectral 
slope distribution from Francis et al. (1992). The arrow marks the number of pairs observed to 
meet the similarity condition. The cases similarity conditions are the same as in Figure 3. 



